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ABSTRACT 

On going searches for supernovae (SNe) at cosmological distances have recently started 
[ to provide large numbers of events with measured redshifts and apparent brightnesses. 

Compared to quasars or galaxies, Type la SNe represent a population of sources with 
q ' well-known intrinsic properties, and could be used to detect gravitational lensing even 

in the absence of multiple or highly distorted images. We investigate the lensing effect of 

■ background SNe due to mass condensations in three popular cold dark matter cosmolo- 
gies (ACDM, OCDM, SCDM), and compute lensing frequencies, rates of SN explosions, 

{Sj . and distributions of arrival time differences and image separations. If dark halos approx- 

imate singular isothermal spheres on galaxy scales and Navarro- Frenk- White profiles on 

O | group/cluster scales, and are distributed in mass according to the Press-Schechter the- 

ory, then about one every 12 SNe at z ~ 1 will be magnified by Am > 0.1 mag (SCDM). 

■ The detection rate of SN la with magnification Am > 0.3 is estimated to be of order 
a few events yr _1 deg -2 at maximum i?-light and Jab < 25, a hundred time smaller 
than the total rate expected at these magnitude levels. In the field, events magnified 
by more than 0.75 mag are 7 times less frequent: about one fifth of them gives rise to 

q \ observable multiple images. While the time delay between the images is shorter than 

H ! 3 days (30 days) in ~ 25% (50%) of the cases (SCDM), a serious bias against the de- 



c/2 



tection of small-separation events in ground-based surveys is caused by the luminosity 
of the foreground lensing galaxy. Because of the flat i^-correction and wide luminosity 
function, Type II SNe dominate the number counts at Jab > 25 and have the largest 
fraction of lensed objects. The optimal survey sensitivity for Type la's magnified by 
Am > 0.75 mag is Jab ~ 23. Magnification bias increases their incidence by a factor of 
50 in samples with Jab < 22, dropping to a factor of 3.5 at 24 mag. At faint magnitudes 
the enhancement is larger for SN II. 
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1. Introduction 

The potential of gravitational lensing as a tool for the determination of cosmological parameters 
or as a probe of cosmogonic models has long been recognized (see Blandford & Narayan 1992 for 
a review). Frequently discussed techniques include detailed models of specific lens systems such as 
multiple quasars (e.g. Falco, Gorenstein, & Shapiro 1991; Grogin & Narayan 1996; Kundic et al. 

1997) , arcs (Lynds & Petrosian 1989; Tyson & Fisher 1995; Tyson, Kochanski, & Dell' Antonio 

1998) , and radio rings (Kochanek et al. 1989), statistical studies of lensed QSOs (Narayan & White 
1988; Fukugita et al. 1992; Kochanek 1996; Wambsganss et al. 1995), and measurements of the 
correlated ellipticities induced by large-scale structure in the images of background galaxies (Babul 
& Lee 1991; Blandford et al. 1991; Kaiser 1992; Villumsen 1996). 

It has been recently realized by many authors (Kolatt & Bartelmann 1998; Metcalf 1999; 
Holz 1998; Marri & Ferrara 1998; Wang 1998; see also Linder, Schneider, & Wagoner 1988) that 
gravitational lensing of distant supernovae (SNe) may provide a new tool for doing cosmology. 
When corrected for light and color curve shapes, Type la SNe can be calibrated to be excellent 
standard candles, with a brightness at maximum l?-light of My = — 19.4±0.15 mag (Riess, Press, &: 
Kirshner 1996; Hamuy et al. 1996). Systematic searches for SNe at high redshifts (Perlmutter et al. 
1998; Garnavich et al. 1998) have already yielded a plethora of events. At this rate of accumulation, 
several hundreds of high-z SNe will be known within the next few years. The uniformity of Type 
la's makes them unique sources for measuring the redshift-luminosity distance relation. 

In the observable (clumpy) universe, gravitational lensing will magnify and demagnify high red- 
shift supernovae relative to the predictions of ideal (homogeneous) reference cosmological models. 
The magnitude and frequency of the effect depend on the power spectrum of density fluctuations, 
the abundance and the clustering properties of virialized clumps, the mass distribution within indi- 
vidual lenses, and the underlying world model. The degree to which weak lensing increases the level 
of noise in the Hubble diagram of Type la SNe, thereby decreasing the precision of cosmological 
parameter determinations, has been discussed by Kantowski, Vaughan, & Branch (1995), Frieman 
(1996), Wambsganss et al. (1997), Holz (1998), and Metcalf (1999). The present paper examines 
instead the high magnification tail of the lensing distribution. The likelihood of a SN at z s ^ 2 
being magnified by Am > 0.3 (> 0.1) mag by a galaxy, group, or cluster is small but non-negligible, 
~ few x 10~ 2 (10%). One can then use recent estimates of the global history of star formation 
to compute the expected frequencies of lensed SNe as a function of cosmic time. The detection 
rate of Type la's with magnification > 0.3 mag is estimated to be of order a few SNe yr _1 deg -2 
at -Tab < 25 mag. Events with Am > 0.75 mag are seven times less frequent (unless specifically 
targeted, Kolatt & Bartelmann 1998). At bright magnitudes, the effect of magnification bias on 
the apparent frequency of lensed SNe is huge. 

On the face of their rarity, which makes a measurement of the magnification probability density 
(hence of the distribution of mass inhomogeneities in the universe) a very difficult task, the modeling 
of individual lensed SN la might have a few advantages over the study of multiple quasars or galaxy 
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correlated ellipticities. Lensed SNe will appear unusually bright for their redshift and will be easily 
identifiable by the proximity of a foreground galaxy, group, or cluster. Their magnification, which 
is related to the integrated mass density along the line-of-sight, is directly measurable to a good 
precision. A supernova that goes off within the Einstein ring of a foreground mass concentration 
may generate multiple images at different positions on the sky. The time delay between the images, 
inversely proportional to the Hubble constant, is longer than a month in about 45 — 55% of the 
cases, and could be measured by comparing the light curves of the two events during a reasonable 
monitoring period. If all the relevant deflection angles to constrain the lens system were known, it 
could be possible to estimate Hq. Small-separation events are associated with galaxy lenses, and 
will remain undetected because of the presence of such a luminous foreground object. 

The plan of the paper is as follows. In § 2 we compute the optical depth for a supernova to 
undergo a lensing event caused by an intervening dark matter halo. We investigate three popular 
hierarchical cosmogonies, two flat models with Qm = 0.35 and 1, and an open model with Q,m = 0.3. 
The population of gravitational lenses is modeled as a collection of singular isothermal spheres on 
galaxy scales and Navarro-Frenk- White halos (Navarro, Frenk, & White 1997; hereafter NFW) on 
group/cluster scales, distributed in mass according to the Press-Schechter theory (Press & Schechter 
1974; hereafter PS). The rates of SN explosions as a function of cosmic time are estimated in § 3, 
together with the expected number counts of magnified events, distributions of time delays and 
image separations, and selection effects. As the rest-frame flux of Type la SNe falls rather steeply 
below 3500 A and we are mostly interested in redshifts z > 1, the counts are computed in the 
/-band. We briefly summarize our results in § 4. 



2. Basic theory 

If the geometry of the universe is well approximated on large scales by the Friedmann- 
Robertson- Walker (FRW) metric, the optical depth for a light-beam emitted by a point source 
at redshift z s and received at z = through a lensing event is (Turner, Ostriker, &; Gott 1984; 
Fukugita et al. 1992) 

t(z s )= dz{l + zf— dMa(M,z,z s )n L (M,z) , (1) 
Jo dz Jo 

where dl/dz = ci7 _1 (l + 2) -1 [fikf(l + z) 3 + 0^(1 + z) 2 + Q^] -1 / 2 is the cosmological line element, 
Qk = 1 — — is the curvature contribution to the present density parameter, a(M, z,z s ) is 
the cross section for lensing as measured on the lens-plane, and i%l(M,z) denotes the comoving 
differential distribution of lenses at redshift z with respect to the variable M (i.e. mass or velocity 
dispersion). Equation (|l]) assumes that each bundle of light rays encounters only one lens, the lens 
population to be randomly distributed, and the resulting r <1. 

The mass distribution in a single lens and the geometry of the source-lens-observer system 
completely determine the lensing cross section. The geometry of the lens system generally af- 
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fects a(M,z,z s ) through the ratio D(0, zi)D(zi, z s )/D(0, z s ), with D(z±,Z2) the angular diameter 
distance between z\ and Z2- In a smooth (homogeneous) FRW universe, one has: 

n / \ c S(x2 -Xi) / n 

D( ^ 2) = ^T^T^) ' (2) 

where S'(x) = sinh(x) if fix > 0, S*(x) = x if = 0, S(x) = sin(x) if £Ik < 0, and 

r ( Z2 dz 
X2 ~ X1 ~ VK J Z1 [n M (i + z)z + n K (i + zy + n A }^ (3) 

with At = 1 if Q,k = and k = \£Ik\ in ah other cases. Here, the density of matter within a cone 
of light rays is always the same as the mean density in the universe (filled beam approach). The 
observed universe, however, is clumpy on galaxy scales, and local mass concentrations are known to 
produce significant changes in the redshift-distance relation (Zel'dovich 1964). Numerical studies 
based on ray-shooting techniques (e.g. Tomita 1998) have shown that a definite distance-redshift 
relation does not exist in a clumpy universe. Analytical approximations are available for light 
beams that subtend a small solid angle, when the probability of intercepting a mass concentration 
is small and light passes through very underdense regions (empty beam approach, Dyer & Roeder 
1973; Ehlers & Schneider 1986). For large apertures (where the mass density in a beam approaches 
the average one), it is widely believed that the angular diameter distance reduces, practically, to the 
FRW case. Note, however, that caustic formation might alter the FRW distance relation even on 
large scales (Ellis, Bassett, & Dunsby 1999, and references therein). For simplicity, in the following 
we will only consider the most widely used filled beam approximation. Depending on the underlying 
cosmology, the discrepancies between the lensing frequences obtained within the various schemes 
for light propagation become significant only for z > 2. 



2.1. Singular isothermal lens 

The standard minimal model to describe the mass distribution of an extended lens (galaxy 
or cluster) is the singular isothermal sphere (SIS), p(r) = a 2 ,/2irGr' 2 , where the mass density p 
as a function of the radial coordinate r is parametrized through a v , the one-dimensional velocity 
dispersion. In the thin lens approximation (e.g. Schneider, Ehlers & Falco 1992; Narayan & 
Bartelmann 1997), this corresponds to a constant deflection angle for incident light rays, (3 = 
An{(j v /c) 2 = 1.4" (0\,/22O km s -1 ) 2 , always pointing towards the lens center of symmetry. For a 
source that, in the absence of lensing, would be seen at an angular distance 9 from this center, the 
lens equation leads to a solution (image) at 0+ = 6 + 6e with magnification p + = 9e/0 + 1, where 
0e = (3Di s /D s is the Einstein radius [hereafter Di s = D(zi, z s ), D s = D{z s ) and D\ = D(zi)]. If the 
alignment is close enough for multiple imaging, 6 < 6e (strong lensing), a second image is produced 
at 9- = 6 — 0e with magnification p^ = Qe/0 — I. The separation between the two images is 29 e, 
and their time delay is cAt = 296 e (1 + zi)DiD s /Di s . At a given angle 9, the total magnification 
of the two images is proportional to their separation. A third image, with zero magnification, 
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is located in correspondence of the lens center, and acquires a finite flux if the singularity that 
characterizes the lens profile is replaced by a core region with finite density. 

A cross section for strong lensing events, csiSj can be easily associated to each SIS lens. 
Measuring cjsis i n angular units, we obtain 

a sls (a v ,z s , Zl )=ne 2 E = 16vr 3 (Eg) . (4) 

The cross section for the image outside the Einstein ring (including both multiply and singly imaged 
objects) being magnified by a factor fi > 1 can be easily related to usis through 

a{n + > n)= ■ ( 5 ) 

Hence the cross section associated with magnification (say) p = 1.316 (corresponding to a bright- 
ening of ~ 0.3 magnitudes) is cr(p+ > 1.316) = IOcjsis- Weakly magnified events are much more 
common, <r(/x+ > 1.1) = lOOusis. In the strong lensing regime, the corresponding value for the 
second (dimmest) image (/i > 0) is 

*{p- > M) = • (6) 

The cross section for a total magnification /i tot = p,+ + p- = 28e/0 > p (with p > 2) is 

o"(Mtot > p) = • (7) 

Therefore, given a population of sources at z s , the distribution of the corresponding magnifications 
(/U + , fi- and /itot) observed at z = is directly proportional to the optical depth t(z s ) when the 
latter is computed by replacing a with Dfasis in equation (|l]). 



2.2. Press-Schechter mass function 

We will follow the approach of Narayan & White (1988) (see also Kochanek 1995b), and use 
the PS theory for our analysis of the abundance of gravitational lenses in a cold dark matter (CDM) 
hierarchical universe. The mass distribution of collapsed halos, assembled through accretion and 
merging processes, follows then directly from the statistical properties of the linear cosmological 
density fluctuation field, assumed to be Gaussian. According to PS, the differential comoving 
number density of dark halos with mass M at redshift z is 



n(M, z) 



1 Po S c (z) 
2nM a 3 



cxp 



M 



2< 



da 2 M 



dM 



(8) 



where po is the present mean density of the universe. The halo abundance is then fully determined 
by the linearly extrapolated (to z = 0) variance of the mass-density field smoothed on the scale M, 
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and by the redshift-dependent critical over density 5 C . The former is determined by assuming 
a post-inflationary spectrum of density fluctuations and by following its linear evolution. In a 
isotropic universe, we can define the power spectrum through the two point correlation function in 
Fourier space (5(k)<5(k')) = (27r) 3 <5z)(k + k')P(fc), where 5d is the Dirac function. The variance is 
then given by the convolution 

1 f°° 

(7 M = ^2j dkk 2 P(k)W 2 (kr ) , (9) 

where W(x) is the top-hat filter function. The critical overdensity is obtained instead by solving 
spherical top-hat collapse and by computing the linear overdensity corresponding to the collapse 
time. Details of these calculations can be found in Lacey & Cole (1993) for vanishing cosmological 
constant, and Eke, Cole & Frenk (1996) for a general flat cosmology. We adopt their analytical 
approach; alternatively one can use iV-body simulations of gravitational clustering to estimate 
5 c (z) as a best fitting parameter (e.g. Lacey & Cole 1994). 

Assuming that every halo virializes to form a (truncated) singular isothermal sphere of velocity 
dispersion a v , mass conservation implies 



a v (M,z) = ^ ro^fA 1/6 



M 



1/6 

(l + z) 1 / 2 , (10) 



with fi(ar) = Q M {1 + zf /[Q, M (l + zf + 9, K {1 + zf + fi A ]. Here r = (3M /4-irp f/ 3 is the 
comoving Lagrangian (i.e. initial, formally corresponding to infinite redshift) radius of the collapsing 
perturbation, z denotes the virialization epoch of the halo, and A (z) is the ratio between its actual 
mean density at virialization and the corresponding critical density for closure of the universe. This 
is determined by assuming that the virialization time is twice the turn-around time (in the limit 
Q,M = 1) one finds the standard result A = 187r 2 ). While in a universe with S7a = the radius 
of a virialized spherical perturbation, r v ; r , is half its turnaround radius, r to , in the presence of a 
cosmological constant the radius at virial equilibrium depends on the vacuum contribution to the 
background energy density, and is always smaller that r to /2 (Lahav et al. 1991). A detailed fitting 
formula for the ratio r Y - lt /r to as a function of and z is given in Kochanek (1995b). 



Equation ( |10] ) relates the PS mass function to the SIS lens profile we discussed in the previous 
section, therefore allowing the computation of the optical depth given in equation (|l|). 



2.3. Navarro-Frenk- White lens 

The number of known gravitationally lensed quasars has grown to around 50, Q and no con- 
firmed lensed system has been detected with angular separation between the images larger than 
6.1". It is well known that models for statistical lensing which combine the PS distribution with 
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SIS profiles tend to overpredict the fraction of large separation lenses (AO ^ 6") with respect to 
the data, even when selection effects are accounted for (e.g. Kochanek 1994; Flores & Primack 
1996). While photometric analyses of E/SO galaxies performed with the Hubble Space Telescope 
(HST) (Tremaine et al. 1994) and a number of gravitational lensing studies (e.g. Wallington & 
Narayan 1993; Kochanek 1995a) indicate that a singular isothermal profile is a good approxima- 
tion for lensing galaxies, high-resolution iV-body simulations have recently shown that virialized 
dark matter halos, formed through hierarchical clustering, have a universal (spherically averaged) 
density profile which is shallower than isothermal (but still diverging like p oc r _1 ) near the halo 
center, and steeper than isothermal (with p oc r -3 ) in its outer regions (NFW), 

P[r) ~ (r/r B )(l + r/r 8 ) 2 ' 1 j 

Here r s is a scale radius, <5nfw is a dimensionless parameter which indicates the characteristic 
density contrast of the halo, and p cr i t is the critical density of the universe^] Measuring the 
halo concentration through the parameter C = r v i r /r s , it follows from the definition of A that 
<5nfw = (^V3)C 3 /[ln(l + C) — C/(l + C)\. The equilibrium density profile of a halo with a given 
virial radius is then completely specified by a single parameter (<5nfw or C)- 

The lens equation for the NFW density distribution has been computed by Bartelmann (1996) 
and Maoz et al. (1997), and can be easily solved adopting standard numerical techniques. The lens- 
ing efficency of a NFW halo can be parameterized by the dimensionless quantity k s = <5NFW/°crit r s/ 
S cr , where S cr = (c 2 /4irG)(D s / D[ s Di) is the critical surface mass density for multiple images. For 
z ~ 0.5, z s ~ 2 and lens masses between 10 12 and 1O 15 M , typical values for k s lie in the interval 
0.07 — 0.4. With respect to SIS profiles containing the same total mass and located at the same 
redshift, NFW lenses have much smaller cross sections for strong lensing (typically ctnfw is a few 
orders of magnitudes smaller than asis , the two cross sections differing less at increasing masses) . 
On the other hand, the average magnification of the images is higher for NFW halos. Thus, de- 
pending on the halo mass (halos of increasing mass being less centrally concentrated, NFW), a 
SIS approximation will overestimate the strong lensing optical depth and underestimate the total 
magnification. 

As suggested by Keeton (1998), a plausible way to reconcile simple statistical lensing models 
with the data on image separations of QSOs is to consider the dissipation and cooling of the 
baryonic protogalactic component and the radial re-distribution of the collisionless dark matter as 
a consequence of baryonic infall (e.g. Blumenthal et al. 1986). The net effect is a compression 
which transforms galaxy sized NFW halos into nearly isothermal distributions. Group and cluster 
halos instead maintain their NFW profile, and are extremely inefficient lenses compared to the 
corresponding SIS. Since baryonic cooling is effective in halos having masses smaller than a threshold 
value M c roughly independent of z (Rees & Ostriker 1977; Silk 1977), the picture adopted in 



5 Note, however, that the slope of the universal density profile in the innermost regions of halos is still subject of 
debate. Moore et al. (1999) find p oc J" -1 ' 5 , while Kravtsov et al. (1998) argue for the presence of a central core. 
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this paper comprises then a collection of isothermal spheres for M < M c and a set of NFW 
halos for M > M c . We have selected the value of M c empirically, by using the Kolmogorov- 
Smirnov test to optimize the agreement with the data on angular separations available in the 
CASTLE (CfA- Arizona-Space- Telescope-LEns) survey. As shown in Figure [l], a 'transition' mass 
of M c ~ 3.5 x 1O 13 M0 provides a best-fit value for all the relevant cosmologies. As multiply imaged 
systems with large angular separations are always produced by massive NFW halos, a reduction 
(with respect to the PS+SIS case discussed in the previous section) of the expected fraction of large 
separation systems to levels compatible with observations is obtained. 

2.4. Lensing optical depth in CDM cosmologies 

We analyze the frequency of lensing events in three different cosmologies with parameters 
suggested by a variety of recent observations: an open model (OCDM, Qm = 0.3, Qa = 0), a 
flat model (ACDM, fl M = 0.35, tt A = 0.65), and an Einstein-de Sitter model (SCDM, Q M = 1, 
= 0). In all cases the amplitude of the power spectrum or, equivalently, the value of the rms mass 
fluctuation in a 8 h^ 1 Mpc sphere, as, has been fixed in order to reproduce the observed abundance 
of rich galaxy clusters in the local universe (e.g. Eke, Cole &: Frenk 1996; Jenkins et al. 1998). [] The 
model parameters are shown in Table |l[ The "cosmic concordance" model proposed by Ostriker 
& Steinhardt (1995) has been selected to represent the whole class of A-dominated cosmogonies 
(ACDM). While this model is able to account for most known observational constraints, ranging 
from globular clusters ages through CMB anisotropics to recent analyses of the Type la SNe Hubble 
diagram, it may be in marginal conflict with the statistics of gravitational lenses (Kochanek 1996). 
Our choice of an open CDM (OCDM) cosmology is guided by the growing evidence in favour of a 
low value of Qm from cluster studies (White & Fabian 1995; Carlberg, Yee, & Ellingson 1997). For 
reference, and as instructive counterpart to the previous variants, we also consider a standard CDM 
(SCDM) scenario. Figure |2| shows the strong lensing optical depth (a = Dfasis for M < M c and 
a = D; 2 ctnfw otherwise) for a source at redshift z s < 3 (also shown, for comparison, are the results 



Here and below we denote the present-day Hubble constant as Ho = lOO/i km s 1 Mpc 



Table 1. Cosmological parameters. 



Model Qm h n as 



SCDM 1.00 0.00 0.50 1.00 0.50 
OCDM 0.30 0.00 0.70 1.00 0.85 
ACDM 0.35 0.65 0.65 0.96 0.87 
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obtained by assuming that every halo relaxes to a SIS). The lensing frequency is the highest in 
SCDM and the lowest (by about a factor of 2 at z s = 3) in OCDM cosmologies, with ACDM models 
in between. In all cases r is dominated by the contribution of SIS lenses. The relative amplitudes 
between different cosmological models depend on a combination of effects: (a) The geometry of 
the lens system affects the optical depth through the ratio (DiDi s /D s ) 2 = (c/Ho) 2 g(zi, z s ) and 
only halos in a narrow redshift interval can be efficient lenses for a source at redshift z s . If z s = 2 
(z s = 1), for example, the maximum contribution to t(z s ) comes from halos at z ~ 0.5 (z ~ 0.3). 
For z\ = 0.5 and z s = 2 (say), the ratios between the geometric factors in the various cosmologies 
are 5scdm/5acdm = 0.58 and ^scdm/^ocdm = 0.85. (b) The path-length to a source at z s = 2 
is (in units of the Hubble radius) 0.71, 0.61, and 0.54 in ACDM, OCDM, and SCDM universes, 
respectively, (c) For a fixed mass, the halo velocity dispersion shows only a weak dependence on 
the world model. At Z\ = 0.5, for example, a^j in ACDM is smaller than in SCDM by 23%. On the 
other hand, OCDM shows the largest velocity dispersions, with larger than in SCDM by 6%. 
(d) The number density of halos obviously influences the lensing frequency. Roughly speaking, for 
objects with mass M <C M*, where M* is such that <tm(-^*) = $c(z), the number of halos contained 
in a Hubble volume is proportional to fijif, i.e. there are approximately three times more lenses in 
SCDM than in low-density universes. By contrast, halos with M > M* form earlier in low-density 
cosmogonies. In SCDM, the formation of the massive (and exponentially rare) objects which are 
responsible for large separation multiple images is delayed to such low redshifts that they fail to 
be efficient lenses for sources at z s > 1. This is graphically shown in Figure ||, where the function 
(M, zi), the mean-free-path per logarithmic mass interval of a beam to encountering 
a lensing event, is plotted as a function of mass for different values of z\ and a source at z s = 2. 
While in the SIS approximation the lensing optical depth is dominated by masses in the range 10 12 
to 10 Mq, i.e. by massive galaxies, groups, and poor clusters, a sharp drop in the mean free path 
is obviously present for M > M c in the combined SIS+NFW case. The most efficient foreground 
lenses are at a redshift less than unity. The larger number of lenses in SCDM more than compensate 
for the longer path lengths to a given redshift and earlier collapse times that characterize open and 
A models. 

It is fair at this stage to point out that, while it has become customary to model (as we have 
done) virialized dark matter halos by smooth distributions, observational data have recently shown 
the importance of considering the presence of baryonic and dark substructure in groups and clusters 
of galaxies (e.g. Williams, Navarro &: Bartelmann 1999; Keeton 1998). The relevance for lensing 
of substructure and, in particular, of massive central galaxies appears to be more cospicuous in 
smaller aggregates. This helps explaining why all the confirmed lens systems detected so far are 
always produced either by galaxies (multiply imaged quasars and radio sources) or by the halos 
of large clusters (arc systems). Those lens systems which are clearly associated with a group or a 
small cluster are always dominated by the lensing effect of a bright galaxy rather than the dark halo 
of the group. While a detailed analysis of the effect of substructures is clearly beyond the scope of 
this paper, we expect that our predictions for groups and small clusters may tend to underestimate 
their lensing optical depth (see also Williams et al. 1999). 
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It is also worth stressing that the PS mass function has been recently shown to overpredict the 
number density of halos in the low-mass tail of the distribution by a factor 1.5-2 with respect to 
high-resolution iV-body simulations (e.g. Gross et al. 1998). In order to test how this inaccuracy 
of the PS formula affects the computation of the strong lensing optical depth, we have re-calculated 
t(z s ) according to equation (|l]) after replacing the PS mass function with the formula derived by 
Sheth & Tormen (1999) in order to reproduce the iV-body outcome. We find that, given the overall 
degree of uncertainty of our calculations, the inaccuracy of the mass function does not represent 
a major issue. For example, in the SCDM model and for z s = 2, the strong lensing optical depth 
derived from iV-body simulations is smaller than its PS counterpart by only 20%. 

Another simplifying assumption we have made involves the spatial distribution of the lens 
population, here modeled with a Poisson process. Indeed, the clustering properties of dark halos 
are known to slightly alter the lensing statistics. In particular, the presence of rich aggregates 
and large voids will increase the frequency of events in the low- and high-magnification tails of the 
distribution (e.g. Jaroszyhski 1992). 



2.5. Galaxy lensing 

As a way to gauge the possible shortcomings of the PS method discussed above [e.g. the 
intrinsic lack of accuracy at the low-mass tail as well as uncertainties in 5 c (z) or in the conversion 
between mass and velocity dispersion or concentration], it is interesting to consider in some detail 
the contribution of known galaxies to the magnification frequency. We model the lensing galaxies 
again as SIS. Their luminosity function, 

= , (12) 



includes just early E/SO types, as late spirals and irregulars are known to make a small contribution 
(< 20%) to the lensing statistics (Kochanek 1996; Maoz & Rix 1993). The velocity dispersion of 
the dark matter in each galaxy halo is related to the luminosity through a Faber-Jackson law 

L = L *{vS ■ (13) 

The assumption that the local properties of E/SO galaxies persist out to z ~ 1 has some observa- 
tional basis (e.g. Lilly et al. 1995; van Dokkum et al. 1998). It also appears that, by and large, 
the statistics of multiply imaged QSOs is consistent with the locally observed galaxy population 
extrapolated to higher redshifts (Kochanek 1993; Maoz & Rix 1993). If the galaxy comoving num- 
ber density does not change with cosmic time, and one neglects all evolutionary corrections to 
equation (13), equation ([[]) can be simplified, and the lensing optical depth becomes proportional 



to a dimensionless constant factor F (Turner, Ostriker, & Gott 1984) 

F = (J^ (j^y J™ dM<r sls (M,z,z s )n L (M,z) = 16vr 3 n* (-J-)' (^) 4 r(2 + a) , (14) 
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where T(x) is Euler's gamma function. We use the best-fit E+SO Schechter parameters from the 
SSRS2 samples (z < 0.05, Marzke et al. 1998): a = -1.00 ±0.09, ra* = (4.4 ±0.8) x 10~ 3 /i 3 Mpc~ 3 , 
and M* = —19.37 ± 0.1 + 5 log h (the absolute blue magnitude associated with a luminosity L*). 
With a fiducial value of er* = 220 kms -1 (Kochanek 1994), this model yields F = 0.017 (cf. 
Fukugita & Turner 1991). For comparison, the all-type luminosity function from the Autofib survey 
(0.02 < z < 0.15, Ellis et al. 1996), when scaled by the Postman & Geller (1984) fraction (31%) of 
E/S0 galaxies, gives a = -1.16±0.05, n* = (7.5±1.3) x 10- 3 h 3 Mpc" 3 , M* = -19.30±0.13+5 log h, 
and F = 0.033. With these assumptions, cosmology only affects the lensing optical depth through 
geometric effects, and the number of lensed sources is the highest in a ACDM model. How does the 
optical depth for strong galaxy lensing compare then with the PS±SIS±NFW frequency shown in 
Figure |2f? For a source at (say) z s = 2, the former (with Ellis et al. parameters) is 1.3 times larger 
than the latter in OCDM, 1.4 times larger in ACDM, and 2.9 times smaller in SCDM. The number 
of galaxy lenses is two times smaller with Marzke et al. parameters. In order to include the effects 
of galaxy groups and clusters, in the following sections we shall estimate the number of magnified 
SNe from the PS formalism in a SCDM cosmology. It is clear from the above discussion, however, 
that the uncertainties in the lensing optical depth are non-negligible. One should also keep in mind 
that the PS mass function, n(M) cx M*> with £ ~ —2, is significantly steeper at the low-mass end 
than that derived from equations (|l2]), and (|l3|), n(M) cx M 2a+1 . Barring selection effects, the 
statistics of gravitational lensing within the PS theory will tend to overestimate the frequency of 
small-separation images. 



3. Supernovae 
3.1. Rates 

The rates of SNe as a function of cosmic time depend on the star formation history of the uni- 
verse, ip(t), the initial mass function (IMF) of stars, (j>(m), and the nature of the binary companion 
in Type la events. The frequency of core-collapse supernovae, SN II and possibly SN Ib/c, which 
have short-lived (^ 50 Myr) progenitors (e.g. Wheeler & Swartz 1993), is essentially proportional 
to the instantaneous stellar birthrate of stars with mass > 8M . For a Salpeter IMF (with lower 
and upper mass cutoffs of 0.1 and 125 M ), the core-collapse supernova rate can be related to the 
star formation rate per unit comoving volume according to 



Is dm <K 



Rn(t) = m^sl Z ^rrT = 0-0074 x 

J 01 amm<p{m) 



M yr- 1 Mpc~ 



yr^Mpc" 3 . (15) 



By contrast, the specific evolutionary history leading to a Type la event is poorly known. SN 
la are believed to result from the explosion of C-0 white dwarfs (WDs) triggered by the accretion 
of material from a companion, possibly another WD or a non-degenerate, evolved star (e.g. Ruiz- 
Lap uente, Canal, & Burkert 1997). The clock for the explosion is set by the lifetime of the primary 
star and, e.g., by the time taken for a non-degenerate companion to evolve and fill its Roche 



-12- 



lobe (or, in the case the companion is another WD, the orbital decay time following gravitational 
wave emission). The evolution of the rate may depend then, among other things, on the unknown 
mass distribution of the secondary binary components and on the distribution of initial orbital 
separations. In this sense Type la SNe follow a slower evolutionary clock than Type IPs, and their 
rate at any given time is sensitive to the past history of star formation in galaxies. Here we follow 
Madau, Delia Valle, & Panagia (1998, hereafter MDP), and parameterize the rate of Type la's in 
terms of a characteristic explosion timescale, t\vd (which defines an explosion probability per WD 
assumed to be independent of time), and an explosion efficiency, rj. The former accounts for the 
time elapsed in the various scenarios from a newly born (primary) WD to the SN explosion itself, 
the latter for the fraction of stars in binary systems that will never undergo a SN la explosion 
because of unfavorable initial conditions. The rate of Type la events at any one time is then given 
by the explosions of all the binary WDs produced in the past that have not gone off already, i.e. 

77 ttdt'Mt') r! dm exp(- f ~ f '~ f ™ U(m) 

RUt) = 1 JJmc r , \ ( 7° — \ (16) 

twd J dmm(p{m) 

where m c = max[3, m(t')], m(t') = (10 Gyr/i') 0,4 is the minimum mass of a star that reaches the 
WD phase at time t', and t m = 10 Gyr/m 2 ' 5 is the standard lifetime of a star of mass m (all stellar 
masses are expressed in solar units). 

Figure || shows the predicted Type la and II rates per unit proper time and unit comoving 
volume for a cosmic star formation history that traces the rise from z = to z ~ 1.5 of the galaxy 
UV emissivity (Madau, Pozzetti, & Dickinson 1998). The stellar evolution model (depicted in 
Figure lb of MDP) assumes a Salpeter IMF and a weakly increasing star formation density above 
z ~ 2, and is able to account for the entire optical background light recorded in the galaxy counts: 
half of the present-day stars - the fraction contained in spheroidal systems - formed at z > 2.5. 
Because of the uncertainties in the amount of starlight that is absorbed by dust and re-radiated 
in the far-IR as a function of cosmic time, these numbers are only meant to be indicative. Our 
estimates may be on the conservative side, as the assumed stellar evolution model produces only a 
fraction ~ 50% of the IR background detected by COBE (Dwek et al. 1998). On the other hand, 
dust along the line-of-sight may affect the detectability of Type II SNe (and also Type la's if the 
explosion delay timescale is short enough that dust is retained in the star- forming regions.) No 
attempt has been made to include the possibility of optically hidden SNe in our calculations. 

The Type la rates plotted in the figure assume characteristic "delay" timescales after the 
collapse of the primary star to a WD equal to tvvd = 0.3, 1, and 3 Gyr, which virtually encompass 
most relevant possibilities. For a fixed initial mass m, the frequency of Type la events peaks at an 
epoch that reflect an "effective" delay At = t\vd + from stellar birth. A prompter (smaller twd ) 
explosion results in a higher SN la rate at early epochs. The explosion efficiency (5 < 77 < 10%) was 
adjusted to reproduce the observed ratio of Type II to Type la rates in the local universe. When 
normalized to the emitted blue luminosity density, the predicted frequencies match rather well the 
data at < z < 0.4 (MDP). 
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3.2. Number counts of SN la 

In this section we will derive the observationally relevant quantity, i.e. the rate of significantly 
magnified SNe in a flux-limited survey, given a cosmological model, lens population, and star 
formation history. An apparent magnitude-redshift relation for unlensed SNe is needed to estimate 
the number counts of distant Type la events. Traditional (i.e. from observed B- to emitted 5-band 
light) if-corrections out to z ~ 0.2 have been calculated by Hamuy et al. (1993) using available 
spectra of nearby SN la. Multi-band photometry has been used to extend these results out to 
z ~ 0.6 (Kim, Goobar & Perlmutter 1996). No accurate ET-corrections have been determined, 
however, for supernovae at z > 0.8, as very few ultraviolet spectra from the local samples with 
good signal-to-noise ratio are currently available. To circumvent this problem we have constructed 
an (unlensed) magnitude-redshift relation for Type la's by using as a template the HST-FOS 
spectrum of SN 1992A taken 5 days after maximum blue light (Kirshner et al. 1993). No correction 
for foreground and intrinsic extinction is apparently needed for this event (Burnstein & Heiles 1984; 
Kirshner et al. 1993). The observed spectrum has been renormalized to a brightness (at maximum 
-B-light) of My = —19.4. The m — z relation in the /-band is shown in Figure | for z < 2.5. Note 
that, because of the steep optical/UV spectrum of 1992A, our template Type la SN at z = 2.5 will 
be much brighter in the near-IR (Jab = 26.9) than in the optical (Iab = 31.2). The X-correction 
significantly weakens the observed / flux relative to the luminosity distance D^ 2 effect. 

For the intrinsic luminosity function of SN la, $>q(I, z), we have adopted a Gaussian distribution 
with standard deviation 0.15 mag (cf. Van den Bergh & McClure 1994). The luminosity function 
of lensed sources, <&(/), can be obtained by convolving &o(I, z) with the magnification probability 
density (Vietri & Ostriker 1983) 

= J dr$o(l',z)V(Am,z) , (17) 

where P(Am, z) is the probability of having magnification Am = I' — I for a source at redshfit z, 
normalized so that P(Am, z)d(Am) = 1. The observed rate on the sky of Type la events is 
then 

where V(z) is the comoving volume element surveyed, and the factor (1 + z) _1 accounts for the 
cosmological time dilation.]] From equation fll7| ) we find 

dt J dZ 1 + Z J-oo J^min(I) 

where log [/i m in(-0] = (I ~ ^Hm)/2.5 is the minimum magnification needed to detect a source in a 
flux-limited survey. The above equation relates the number counts to the probability distribution 



7 Note that, as in a snapshot survey SNe are not necessarily detected at maximum light, the total number of events 
will actually depend on the effective timescale over which a SN can be observed above threshold. 
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of magnification discussed in § 2 for SIS lenses. In the case of NFW halos, we first solve the 
lens equation for a set of values of k s and of the dimensionless impact parameter y = D\Q jr s . The 
magnitudes of the resulting images and the values of the corresponding impact parameters are stored 
in a file. The optical depth for a given magnification and source redshift is then computed using 
equation ([!]): the suitable cross sections are calculated by interpolating over the list of solutions 
of the lens equation, and by using the method described in the appendix of NFW to evaluate the 
scale radius r s of a selected halo. 

To compute the total number of supernovae brighter than Ii; m that have been gravitationally 
magnified by more than Am magnitudes, the function V(/i) in equation (|l9| ) must be integrated 
from max (/u m i n , 10 Am//2 ' 5 ) to infinity. Our treatment will be limited to magnifications Am > 0.3 
mag, i.e. to images that are more than 2a away from the mean intrinsic luminosity of SN la. In this 
regime, the weak (with typical magnifications of the order of a few per cent) lensing due to large 
scale structure can safely be neglected. Figure |6| shows the counts of SNe la at maximum B-light 
for a SCDM cosmology. The rate of detection in the absence of gravitational lensing effects (which 
is very close to the total expected number of events, since the probability of lensing is small) is 
compared with the estimated rates of multiply-imaged events (obtained by requiring that at least 
two images are sufficiently bright to enter the flux-limited sample; e.g., for SIS lenses, using eq. |6|), 
and with the rates of individual SNe magnified by more than 0.3 and 0.75 magnitudes (for strongly 
lensed SNe, only the brightest image is considered; i.e., for SIS lenses, V(fJ>) is computed using the 
cross section in eq. || while, for NFW halos, the above mentioned numerical scheme is adopted). 
Note that, for low magnifications (/U ^ 2), the lensing cross section of NFW halos is comparable 
(and generally larger) than that of the corresponding SIS halos. Only considering multiply imaged 
systems, the relation cjnfw "C cxsis holds. 

In a one square degree field, and for an effective observation duration of 1 yr, the total expected 
number of Type la SNe down to Jab < 25 mag is (assuming twd = 1 Gyr) ~ 390. About 3.7 events 
will be gravitationally lensed by more than 0.3 mag, 0.5 by more than 0.75 mag, and only 0.1 will 
have detectable double images. These numbers are sensitive to the assumed explosion timescale. 
Given a well-determined cosmic star formation history, it may be possible to infer something about 
the nature of Type la SNe from the observed number counts. It is also interesting to note that, 
in general, gravitational magnification facilitates the detection of SNe at higher redshift in a flux- 
limited sample, and the fraction of lensed SNe tends to be higher at fainter magnitudes if the rate 
of events increases with look-back time and the gain in available volume is not balanced by a large 
ET-correction. For Am > 0.75, this is true up to Jab ~ 23. For fainter magnitudes, the increase 
in physical volume is smaller and the SN rate saturates (cf. Fig. |I|); deeper searches will not lead 
then to an increased fraction of magnified sources. This defines an optimal survey sensitivity in 
the I-band for Type la's magnified by Am > 0.75 mag of Iab ~ 23 (cf. Linder et al. 1988). By 
contrast, there is no optimal flux threshold for detecting double images (Fig. |6|d) in succession. 
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3.3. Number counts of SNe II 

Type II SNe are a couple of magnitudes fainter than SN la in the rest-frame blue band. 
Contrary to Type la's, they have a very wide luminosity function, which is well represented by a 
broad Gaussian distribution with (Ms(max)) = —17.2 (h = 0.5) and dispersion 1.2 mag (Tammann 
& Schroder 1990).^] The spectra of most SN II at maximum light approximate well a single- 
temperature Planck function of about 15000 K from the UV through IR - by contrast, Type la's 
exhibit a prominent UV deficit relative to the blackbody fit at optical wavelengths (e.g. Filippenko 
1997). The difference in behavior of the i^-correction is clearly visible in Figure [?], where the /-band 
m — z relation for Type II's is plotted for three different photospheric temperatures at peak. In 
the absence of dust reddening, the intrinsically fainter SN II appear brighter than SN la already at 
z ~ 1.5. This has a strong effect on the expected number counts, shown in Figure || At Jab < 24 
mag, there is one SN II every 2.3 Type la's. The detection rate increases to one Type II every 
1.2 SN la at Jab < 25, and to one every 0.56 at Iab < 26. Because of the flat iT-correction 
and wide luminosity function, the fraction of lensed events is higher for Type II's, as only small 
magnifications are needed to bring high redshift events above the flux threshold (see also Linder et 
al. 1988). 



3.4. Multiple images and time delays 

Although exceedingly rare in the field, multiply imaged Type la SNe are of particular interest 
for cosmography, as the propagation time from the source to the observer varies from one image to 
another, and can be accurately measured if the source has a well defined light curve. The time delay 
At is proportional to the absolute scale of the lens geometry and hence inversely proportional to 
the Hubble constant (Refsdal 1964). Given a reliable time delay measurement, a well-constrained 
lensing mass distribution (with a measured redshift and velocity dispersion), and a world model, 
one could estimate Hq. 

The cumulative probability functions of image separations AO and time delays At for a SN 
at z s = 2 in the field are shown in Figure (||). About 60% (SCDM) of double images will have 
angular separations < 1", and will therefore be flagged as multiple events in ground-based surveys 
only if the time delay is comparable or larger than the SN decay time. The fraction of images with 
A9 > 4" is of order 1.5%. Barring selection effects (see below), approximately 25% of mirror events 
will be observed with an arrival time difference of < 3 days, 40% with < 10 days, 50% with < 30 
days. The frequency distributions have been computed by modifying equation ([[]) to account for 
that fact that, at any redshift, only halos having mass larger than a threshold value can generate 
two images with a given angular separation or time delay. Normalizing the total probability to 



8 For the present purposes, we shall ignore the possibility of an extended tail of Type II events at bright (Patat et 
al. 1994) and faint (Woltjer 1997) magnitudes. 
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unity, and considering only SIS halos with M < M c (whose contribution dominates the strong 
lensing optical depth), one gets 



V{> A9\z s ) =t- 1 {z s ) / dz(l + z) 



dl 



Mc 



dMDfasisn L , 



(20) 



dz JM(Ae,z,z s ) 

where M is the mass that, for given z and z s , produces a separation AB [equal to 87r(a v (M)/c) 2 
(Dis/D s ) for a SIS]. Similarly, one obtains 



V(> At\z s ) = T -\z,) / dz(l + zf — 



dl 
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Dfasis n L , (21) 



where At max is the time delay at B = Be and M solves At max (M, z,z s ) = At. Note that, because 
of the small strong lensing cross section of NFW halos, the presence of very large dark matter halos 
in OCDM and ACDM universes does not create an extended tail towards large separations and 
long delays. 



3.5. Selection effects 



A peculiar feature of SNe la relative to QSOs is the relatively short lifetime of their emission. 
While, for example, the finite angular resolution 9 rcs of the observations may not allow to distinguish 
multiple images with separation on the sky A9 < 9 res , if the time delay between different images is 
comparable or longer than the typical decay time of a SN light curve [£sn ~ 20(1 + z s ) days], two 
mirror copies will appear roughly in the same position of the sky but at different times. In this 
sense, even two images with AO < 9 rcs could in practice flag a strong lensing event if At > tsw. the 
same SN would go off in two (or more) incarnations of the same explosion, and the different images 
will be seen in succession. Because of the geometry of the source-lens-observer system, longer time 
delays will be obtained when the deflector is closer (in redshift) to the radiation source. Consider, 
e.g, a SN at z s = 2 in a SCDM cosmology, lensed by a massive SIS halo (M = 1.2 x 1O 13 M , 
<r v = 250 kms" 1 ) at z t = 1.2. Two images will result, with AB = 0.83" < 9 TCS ~ 1" (say). For 
a source-lens alignment 9 = Be/2, the magnifications of the images will be /x+ = 3 and ji- = 1, 
and the arrival time difference At = 81 days. Unresolved images with shorter delays may still 
be recognizable due to the large discrepancy between the apparent expected magnitude and the 
one observed, and a characteristic doubly-humped light curve. At fixed z\ and z s , the angular 
separation between the lensed images depends only on the mass of the deflector, while the time 
delay is anti-correlated with the magnification: smaller delays correspond to more perfect source- 
lens alignments and thus to larger magnifications. Nearly synchronous images will then have rather 
similar brightnesses. A SN at z s = 2, lensed by a SIS halo at z\ = 0.5 with mass 8.2 x 10 12 M@ 
(a v = 220 kms -1 , SCDM), will produce when B = 0.15 Be two well resolved images (AB = 1.6"), 
with magnifications /x + = 7.8 and /x_ = 5.8. The arrival time difference is only 20 days in this 
example. For SIS lenses, the maximum time delay occurs when B = Be (At max oc 9 E ), hence less 
massive deflectors will always produce small-separation events and relatively short delays. 
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A serious bias against the detection of small-separation images is due to the luminosity of the 
lensing galaxy. Galaxy-type halos with velocity dispersion a v < 200 kms -1 completely dominate 
the rate of small-separation events. These must then be detected in the background of a lens 



elliptical galaxy, most probably at z\ ~ 0.5. Figure 10 compares the observed surface brightness of 
a foreground -L* elliptical with a background lensed Type la SN as a function of angular distance 
from the galaxy center. The effect of atmospheric seeing has been approximated by convolving 
the true source flux with a Gaussian of dispersion a = FWHM / y/8 In 2. Even in the case of 
excellent FWHM= 0.6" seeing, the lens will greatly outshine the faintest SN image. Double events 
will only be detectable for nearly perfect lens-source alignments, an extremely rare case. In more 
modest seeing conditions (FWHM= 1"), even the brightest image may remain undetected. From 
the ground, multiple SN events will be more easily observed in the background of a galaxy group 
or cluster. 

In a magnitude-limited sample, the sources that are observed to be gravitationally lensed 
include not only the lensed objects that are intrinsically brighter than the flux limit, but also sources 
that are intrinsically fainter and are brought into the sample by the magnification. Relative to the 
lensing optical depths shown in Figure lensed systems will then be overrepresented among SNe of 
a given apparent magnitude. We have quantified this effect from the Type la and Type II number 
counts, by defining the magnification bias B(< /n m ) as the ratio between the actual flux- limited 
counts of sources magnified by more than Am and the counts of lensed SNe that are intrinsically 
brighter than the flux limit, 

dV(z) R(z 



B{< hm) - , v( v R( v h , (22) 

J dz 1 + zJ-oc UV ' A0Am/2.5 P ' 

where = max[// m i n (/), 10 Am//2 5 ]. This definition generalizes that given in Turner, Ostriker, & 
Gott (1984) and Fukugita & Turner (1991) by taking into account the redshift dependence of the 
lensing optical depthj^j The magnification bias corresponding to Am = 0.3 and 0.75 mag is plotted 
in Figure for Type la (r WD = 1 Gyr) and Type II SNe (T = 15000 K). At J AB < 22 the effect is 
huge, as the apparent frequency of events with Am > 0.75 is boosted by a factor of 35 — 50. Note 
that at these magnitude levels the counts are still nearly Euclidean: the large bias is actually due 
to the rapidly increasing lensing frequency with lookback time. In the standard treatment with 
non-evolving optical depth, a large magnification bias is typically associated with integral number 
counts that increase faster than f~ 2 , where / is the observed flux (e.g. Schneider, Ehlers & Falco 
1992). 



In the notation of Fukugita & Turner, our approach is equivalent to a flux-dependent lensing optical depth, as 
fainter objects are typically at higher redshifts than bright ones. 
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3.6. Microlensing 

A modulation of the SN light curve may take place if, on top of the galaxy, group, or cluster 
magnification, a microlensing event occurs due to individual stars or MACHOs in the galaxy halo 
or intracluster medium. The timescale for microlensing-induced variations is determined by two 
competing effects: the expansion velocity of the SN photosphere and the relative transverse motion 
between the microlens and the light bundle. The radius of a SN at maximum light is r§N ~ 10 15 
cm, corresponding to an observed angular dimension of #sn = ^sn/-Cs- A compact object along the 
line-of-sight to the SN will give rise to a detectable time dependent amplification only when 9e^> 
#sn ~ |#| (Schneider & Wagoner 1987), where is the angular impact parameter. Modeling the 
deflector with a Schwarzschild lens of mass M (i.e. neglecting the shear induced by the gravitational 
field of the host halo) , the first condition gives 

M » M E = 1.83 x 1(T 4 M h ( —?-±- \ . (23) 

\H D s DisJ 

Less massive MACHOs can efficiently magnify the SN only during its fast light-rising time. For 
M 3> Me, the duration of the microlensing event roughly corresponds to the time needed by the 
photosphere to fill the Einstein ring of the lens (see also Gould 1994). For a typical photospheric 
expansion velocity of f» 10 9 cm s" 1 , one obtains 

This is usually much shorter than the characteristic crossing time of the Einstein ring 
* u-i/2 { HoDlDls \ 1/2 { M NVV gOOkms-^ \ 

tEK3Jyih 1 {-^DT) U^MoJ { v ± ) {1 + Zlh (25) 

where v± ~ \pio v is the relative transverse velocity between microlens and source. In a SCDM 
cosmology with z s = 2 and z\ = 0.5, one finds i p h ~ 49.2 days (M/1O~ 3 M ) 1 / 2 and tE ~ 
2.9 yr (M/1O~ 3 M ) 1//2 (2OO kms^ 1 /v). The microlensing timescale is then comparable with i$N ~ 
20(l+z s ) days (SN la) for very small microlens masses. Only MACHOs with 1(T 4 < M/M Q < 10~ 3 
can then produce detectable 'noise' that will degrade the light curves of the macroimages. More mas- 
sive objects will boost the magnification of the SN image for the entire duration of the macroevent. 



4. Summary 

All previous studies of the influence that a population of deflectors has on the observational 
properties of background objects have considered sources like galaxies or quasars that are quite 
inhomogeneous, i.e. have a wide range of intrinsic luminosities, ellipticities, etc. By contrast, 
Type la SNe represent a population of objects with well-known intrinsic properties, and could - 
in addition to their use as reliable cosmological distance indicators and valuable tracers of the 
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star formation history in galaxies - potentially flag the occurrence of a gravitational lensing event 
even in the absence of multiple or highly distorted images. We have investigated the lensing 
effect of background SNe due to virialized dark matter halos in CDM cosmogonies, and computed 
lensing frequencies, rates of SN explosions, and distributions of arrival time differences and image 
separations. Within the assumption that dark halos approximate singular isothermal spheres on 
galaxy scales and NFW profiles on group/cluster scales, and are distributed in mass according to 
the Press-Schechter theory, about one every 12 Type la SNe at z ~ 1 will be magnified by Am > 0.1 
(SCDM). Using recent estimates of the global history of star formation to compute the expected 
frequencies of SNe as a function of cosmic time, we have derived a detection rate of Type la's with 
magnification Am > 0.3 mag of a few events yr _1 deg -2 at maximum .B-light and Jab < 25. Events 
with Am > 0.75 mag are 7 times less frequent. Only about one fifth of them give rise to observable 
multiple images: the SN will appear as multiple events with identical light curves, separated in 
time and differing only by the scaling of their amplitudes. Because of the flat if-correction and 
wide luminosity function, we find that Type II SNe will dominate the number counts at Jab > 25 
and have the largest fraction of lensed objects. At bright magnitudes, the effect of magnification 
bias on the apparent frequency of lensed SNe with Am > 0.75 is huge, a factor of 35 — 50 in samples 
with Jab < 22. The enhancement of the lensing probability beyond the optical depth expectation 
drops to a factor of 3.5 at 24 mag (26 mag) in the case of Type la SNe (SN II). The apparent 
magnitude of the lens galaxy introduces a serious selection effect, as it reduces the detectability of 
small-separation multiple events in ground-based searches. 

On the face of their rarity, it is not unconceivable that large consortia of new and existing 
SN search teams may conduct in the near future deep pencil beam surveys and be able to detect 
magnified Type la's up to z ~ 2 (Wang 1998). Current searches have been limited so far to R ~ 24 
mag. A second-epoch deep (Jab ~ 27 mag) observation of the Hubble Deep Field by Gilliland, 
Nugent, & Phillips (1999) has recently revealed a z ~ 1.3 Type la SN. The combination of a new 
generation of large telescopes with 1 degree field of view and multi-fiber spectrometers, together 
with faster techniques such as the snapshot distance method (Riess et al. 1998), could soon yield 
more than ~ 500 new events per year. As suggested by Kolatt & Bartelmann (1998), one could also 
look for lensed SNe in the background of massive clusters. A nearby rich cluster at z\ = 0.1 with 
virial mass of 10 15 Mq and NFW profile will magnify by more than 0.3 mag every point source lying 
at less than 93" from its center. The expected rate of SNe la behind such a cluster with a mismatch 
of > 0.3 mag is then ~ 1.2 events per year with Jab < 26 (SCDM, twd = 1 Gyr). Along the same 
lines, Kovner & Paczynski (1988) have proposed the possible identification of Type II SNe in the 
lensed giant arcs observed in the cores of clusters of galaxies. If a SN goes off within the caustic and 
close to the cusp of the cluster potential, it will be detected as three (or more) spatially-separated 
events with time delays as short as a few weeks or even days. The blue color of the arcs suggests 
significant star formation hence a high rate of SN explosions. Incidentally, if the cluster potential 
were understood well enough that a reasonable estimate of the time delay was available, one could 
exploit the occurrence of delayed mirror events to study such transient phenomena in much greater 
details that would be possible without prior warning. 
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Looking further into the future, it is interesting to consider the detectability of lensed Type 
la and Type II SNe with the Next Generation Space Telescope (NGST). A typical SN II at z = 5 
would give rise to an observed flux of 8 nJy (SCDM) at 2.6 /j,m. At this wavelength, the imaging 
sensitivity of an 8m NGST is around 1.5 nJy (10 4 s exposure and 10a detection threshold), while the 
moderate resolution (A/AA = 1000) spectroscopic limit is about 50 times higher (10 5 s exposure per 
resolution element and IOct detection threshold, Stockman et al. 1998). Depending on the history 
of star formation at high redshifts, the NGST could detect ~ 10 Type II SNe (at peak brightness) 
per 4' x 4' field per year at z > 4 (MDP; cf. Dahlen & Fransson 1999; Miralda-Escude & Rees 
1997). (In principle, Type la's could also be imaged at extreme redshifts, z > 10, but the long delay 
times from stellar birth characteristic of these events make a detection rather improbable.) We find 
that, even for sources at these early epochs, the optical depth for strong lensing is small, r ~ 0.008 
for z s = 7 and SCDM. Marri & Ferrara (1998) have recently explored in details the observational 
perspectives for Pop III SNe detection with the NGST when gravitational magnification is taken into 
account. They computed lensing frequencies and magnification maps using a ray-shooting numerical 
technique and assuming point-like lenses. Our more realistic treatment implies probabilities that 
are much lower than theirs, as isothermal spheres and NFW halos are very inefficient magnifiers 
compared to point-like deflectors of the same mass. 

We have benefited from useful discussions with B. Bassett, R. Blandford, C. Keeton, N. Pana- 
gia, and P. Schneider. We thank R. Kirshner for providing us with the HST-FOS spectrum of SN 
1992A, and the referee, C. Kochanek, for many comments that vastly improved the manuscript. 
Support for this work was provided by NASA through ATP grant NAG5-4236. 
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Fig. 1. — Predicted probability distribution of angular separations between the two brightest images 
of a lensed point source at z s = 2 in different cosmologies. The lenses are approximated as singular 
isothermal spheres for M < 3.5 x 1O 13 M , as Navarro-Frenk- White halos otherwise, and are 
distributed in mass according to the Press-Schechter theory. The angular selection effects have 
been included following Kochanek (1999, private communication). The data with errors bars are 
taken from the CASTLE survey. 
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Fig. 2. — Strong lensing optical depths for a point source at z s in three different hierarchical 
cosmogonies. The mass distribution of the lenses is described by the Press-Schechter theory. Solid 
lines: a lens having mass M is modeled by a singular isothermal sphere for M < 3.5 x 1O 13 M , and 
by a NFW profile otherwise. Dashed lines: each lens is modeled by a singular isothermal sphere. 
In this case the lensing optical depth for magnifications Am > 0.3 (> 0.1) mag is 10 (100) times 
larger than plotted here. 
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Fig. 3. — Mean-free-path per logarithmic mass interval of a beam to encountering a strong lensing 
event, as a function of mass for different values of the lens redshift zi and for different cosmological 
models. The source is assumed to be at z s = 2. Each lens is modeled by a singular isothermal 
sphere. 
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Fig. 4. — Estimated Type la and II (rest-frame) frequencies as a function of redshift. Dot-dashed 
line: SN II rate. Dotted line: SN la rate with twd = 0.3 Gyr. Solid line: SN la rate with t\vd = 1 
Gyr. Dashed line: SN la rate with twd = 3 Gyr. The model assumes a SCDM cosmology. The 
data points with error bars have been derived from the measurements of Cappellaro et al. (1997), 
Tammann et al. (1994), Evans et al. (1989), and Pain et al. (1997). See MDP for details. 
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Fig. 5. — /-band magnitude-redshift relation (SCDM) for Type la SNe. The apparent AB mag- 
nitude has been computed by redshifting and dimming the rest-frame H ST-FOS spectrum of SN 
1992A, and renormalizing it to My = —19.4. (The renormalized spectrum is shown in the inset, 
see Kirshner et al. 1993). The emitted power L v is in ergss -1 Hz -1 . 




Fig. 6. — Differential number counts of Type la SNe at maximum 5-light versus limiting /-band 
magnitude of the survey (SCDM universe), a) SNe in the absence of lensing; b) lensed SNe with 
magnification greater than 0.3 magnitudes (for strongly lensed sources only the brightest image is 
considered); c) lensed SNe with magnification greater than 0.75 mag (for strongly lensed sources 
only the brightest image is considered) ; d) strongly lensed SNe with two images above the detection 
threshold (only the dimmest of the detectable images is considered). Dotted line: rate assumes a 
time delay twd = 0.3 Gyr. Solid line: twd = 1 Gyr. Dashed line: t\vd = 3 Gyr. The effect of 
dust extinction on the detectability of Type la's has not been included in the models. 
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Fig. 7. — Magnitude-redshift relation (SCDM) for Type II SNe. The apparent / magnitude has 
been computed in the absence of dust reddening, and by assuming a blackbody spectrum at different 
temperatures. In all cases the emitted power has been renormalized to Mb = —17.2. Dashed line: 
T = 10000 K. Solid line: T = 15000 K. Dotted line: T = 20000 K. 
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Fig. 8. — Differential number counts of Type II SNe at maximum l?-light versus limiting /-band 
magnitude of the survey (SCDM universe, z max = 5). a) SNe in the absence of lensing; b) lensed 
SNe with magnification greater than 0.3 magnitudes (for strongly lensed sources only the brightest 
image is considered); c) lensed SNe with magnification greater than 0.75 mag (for strongly lensed 
sources only the brightest image is considered); d) strongly lensed SNe with two images above 
the detection threshold (only the dimmest of the detectable images is considered). Dashed line: 
rate assumes a photospheric temperature of T = 10000 K. Solid line: T = 15000 K. Dotted line: 
T = 20000 K. 
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Fig. 9. — Normalized cumulative probability distributions of angular separations and time delays 
between the two brightest images of an SN at z s = 2. Solid line: SCDM. Long-dashed line: ACDM. 
Short-dashed line: OCDM. 
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Fig. 10. — Observed /-band light profile for an elliptical galaxy at z\ = 0.5 (with typical luminosity 
and color, Mb = —19.37 + 5 log h, rest-frame B — V = 0.95), as a function of angular distance from 
the galaxy center (solid line). A de Vaucouleurs' law has been adopted with effective radius R e = 4 
kpc. The galaxy acts as a SIS lens with o~ v = 225 kms" 1 for a background Type la SN at z s = 1.5 
(0e = 0.83"). Dashed lines: peak surface brightness of the dimmest (0- = 9 — 9e, left curve) 
and brightest (#+ = 9 + 9e, right curve) SN images for different source-lens alignments 6. The 
filled points along the curves identify multiple events corresponding to the same impact parameter. 
Only one image is produced at 0+ > 29 e- The empty square corresponds to a single image with 
Am = 0.3 mag. The light profiles of SN and galaxy have been convolved with a Gaussian point 
spread function of FWHM 0.6". Dotted line: SN brightness in the absence of lensing. A SCDM 
cosmology has been assumed. 
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Fig. 11. — Magnification bias versus apparent limiting magnitude. Solid lines: Type la SNe with 
t\vd = 1 Gyr. Dashed lines: Type II SNe with T = 15000 K. The two set of curves have been 
computed for magnifications Am > 0.3 and Am > 0.75 mag. For strongly lensed sources, only the 
brightest image is considered. 



